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Abstract:  In  this  paper  we  apply  vector  s-algorithm  to  solve  a  system  of  linear  equations 
arising  in  the  method  of  moments  solution  of  electromagnetic  scattering  from  wire 
antennas.  The  method  does  not  require  the  inversion  of  a  matrix.  Hence,  it  avoids  the 
problems  associated  with  matrix  inversion  of  singular  or  ill-conditioned  matrices. 
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1.  Introduction 

This  paper  presents  an  iterative  technique  for  solving  the  system  of  equations  obtained  in  a 
method  of  moments  solution.  The  method  is  based  on  vector  s-algorithm  (VEA)  [1].  The 
application  of  the  method  is  illustrated  by  applying  it  to  find  the  current  induced  on  a  wire 
excited  by  a  normally  incident  plane  wave.  Matrix  inversion  is  not  required  as  the  method 
is  iterative.  The  solution  process  can  be  stopped  when  a  pre-defined  convergence  criterion 
is  met.  This  allows  the  user  to  specify  the  degree  of  accuracy  in  the  solution.  Some  of  the 
iterative  techniques  utilized  so  far  include  the  conjugate  gradient  method  [2]  and  the 
method  of  steepest  descent  [3].  VEA  presents  an  alternative  to  solving  a  system  of 
equations,  especially  when  other  methods  fail  due  to  ill-conditioned  or  a  singular  matrix 
[4], 


2.  Vector  s  -  Algorithm  (VEA) 
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Consider  the  solution  of  the  following  system  of  linear  equations: 

AX=b  (1) 

where  A  is  a  known  NxN  moment  method  matrix,  A  =  \ari  I ;  X  is  the  unknown  Nx  1 
column  vector, X  =  (xl,x2,...xN)T ,  and  b  is  known  Nx  1  forcing  function  vector, 

b  =  {bl,b2  ...bN)  .  In  VEA,  estimates  of  X,  designated  X  ,  are  obtained  without  the  use  of 
matrix  inversion. 

The  starting  vectors  }  and  {s1^’ },  each  of  length  N,  are  constructed  from  the  initial 
values: 

{s“)={0),  q  =  1,2,...  (2) 

hrMs|0|Mo} 
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The  vectors,  S ^  •••  sn^)  >  can  be  calculated  using  Gauss-Seidel  relaxation  [5]: 
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i  =  1,  2 ,...N 

q  =  1, 2, ... 


(4) 


Once  we  have  at  least  two  initial  vectors  S 1 0 1  and  S 1 11 ,  then  the  VEA  can  be  started  to 
obtain  (r  +  V)th  order  iterates,  which  are  given  by: 


q  =  0, 1,  2,  . . . 

f  =  l,  2, ...  V 


(5) 


The  value  of  index  r  in  the  above  equation  takes  on  values  of  0,  1,  2,  ....  In  (5),  the  even 
order  vectors  •••  provide  estimate  of  the  solution  vector  A.  The  computations 

are  stopped  when  the  following  condition  is  met: 

-b) 

— ^  <  NCf  (6) 

2 


where  C  f  is  a  pre-defined  convergence  factor. 


3.  Numerical  Results 
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We  show  the  application  of  VEA  in  solving  the  current  induced  on  a  thin  wire  scatterer. 
The  integral  equation  for  the  induced  current,  given  by  Pocklington’s  equation,  is  solved 
using  pulse  expansion  and  triangular  testing  functions  [6] . 

Fig.l  shows  the  current  induced  by  a  normally  incident  plane  wave  ( E0  =  1  V/m)  on  a  wire 
scatterer  of  length  £  =  ^ ,  radius  a  =  /^qqq  >  frequency  /  =  300  MHz  and  convergence 

factor  Cf  =  0.01.  The  number  of  pulse  expansion  functions,  A,  is  31.  We  have  run  the  code 

for  higher  values  of  N  and  obtained  results  in  agreement  with  those  obtained  from  direct 
methods.  Presently  work  is  in  progress  to  apply  the  algorithm  to  solve  2-dimensional  (flat 
plate)  and  3-dimensional  (conducting  sphere)  scattering  problems  using  triangular  patches 
in  the  method  of  moments.. 


4.  Conclusion 

In  this  work,  we  have  shown  that  vector  s-algorithm  can  be  utilized  to  solve  a  system  of 
linear  equations  obtained  in  a  moment  method  solution  of  scattering  problems.  The  results 
for  the  computation  of  current  distribution  on  a  thin  wire  scatterer  show  very  good 
agreement  with  direct  method  (matrix  inversion).  Since  this  method  is  iterative,  it  does  not 
require  inversion  of  a  matrix.  This  can  be  advantageous  if  one  encounters  a  singular  or  ill- 
conditioned  moment  method  matrix. 
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Figure  1:  The  current  induced  on  a  thin  wire  of  length  =  X/2  ,  radius  =  )J  1 000  and  N  =  31. 
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